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It has been predicted theoretically and observed experimentally that disorder leads to spatial fluctu- 
ations in the superconducting (SC) gap. Areas where SC correlations are finite, coined SC islands, 
were shown experimentally to persist into the insulating side of the superconductor-insulator tran- 
sition. The existence of such (possibly weakly coupled) SC islands in amorphous thin films of 
superconducting material accounts for numerous experimental findings related to superconductor- 
insulator transition and non-monotonic magneto-resistance behavior in the insulating region. In this 
work, a detailed analysis pertaining to the occurrence of SC islands in disordered two-dimensional 
superconductors is presented. Using a locally self-consistent numerical solution of the Bogoliubov- 
de-Gennes equations, the formation of SC islands is demonstrated, and their evolution with an 
applied perpendicular magnetic field is studied in some detail, along with the disorder-induced 
vortex-pinning. While mean-field theory cannot, in principle, explore phase correlations between 
different islands, it is demonstrated that by inspecting the effect of a parallel magnetic field, one 
can show that the islands are indeed uncorrelated SC domains. Experimental predictions based on 
this analysis are presented. 

PACS numbers: 74.20.-z, 74.45. +c,74.81.-g 



I. INTRODUCTION 

Interplay between disorder and superconductivity has 
been at the focus of attention for quite a long time. 
More than four decades ago, it has been established 
that the effect of weak disorder on superconductiv- 
ity is not substantial. 1,2 Strong disorder, however, may 
have a profound effect, driving the system from a su- 
perconducting (SC) to an insulating state. Such a 
superconductor-insulator transition (SIT) was observed 
in two-dimensional amorphous SC films. 3 Upon decreas- 
ing the film thickness or increasing a perpendicular mag- 
netic field, these films (which are held below their bulk 
critical temperature) exhibit a transition from a SC state, 
characterized by a vanishing resistance as T — > T c , to 
an insulating state, in which the resistance diverges as 
T — > 0. The possibility of tuning the system continu- 
ously between these two extreme phases may be a man- 
ifestation of a quantum phase transition (which, strictly 
speaking, occurs at zero temperature). 

Despite the substantial amount of experimental data, 
and numerous theoretical investigations, there are still 
several unresolved issues pertaining to the physics of 
these systems. One of the main puzzles is the mechanism 
by which the magnetic field destroys the SC correlations. 
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The "dirty boson" theory, 5,6 describing bosons (Cooper 
pairs) in a random potential, regards the SIT as a transi- 
tion into a Bosc-glass phase, in which the pairing ampli- 
tude is finite in the sample but its phase is strongly fluc- 
tuating, giving rise to an insulating state. Such a phase 
is characterized by a pair- vortex duality, and its hallmark 
is a universal value for the resistance Rq = h/ie 2 at the 
transition. While some experiments are consistent with 
that prediction, 7-9 others are not. 10 ' 11 Beside the issue of 
universal resistance, there are other observations which 
cannot be simply explained by the dirty boson model, 
for instance, the temperature dependence of the crossing 
point and the classical XY critical exponent. 12 

The role of the large fluctuations in the local SC gap, 
A, and the formation of SC islands have been empha- 
sized in several other works. 13-16 In Rcf. 13 the authors 
predict the existence of such islands by calculating the 
mesoscopic fluctuations of the order-parameter at the 
mean-field level. Since phase fluctuations are ignored by 
mean- field calculations, it was pointed out that there is 
no real SC transition within that approximation. Later 
it has been shown 14 that inclusion of quantum fluctua- 
tions (beyond mean field) indeed results in a quantum 
phase-transition. A detailed mean-field study of SC cor- 
relations in disordered two-dimensional films has been 
presented in Ref. 15, employing a locally self-consistent 
solution of the mean-field Bogoliubov-de Gennes (BdG) 
equations 17 at finite temperature and at zero magnetic 
field. It was indeed found that disorder induces strong 
fluctuations of A in space, resulting in regions with large 
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order parameter (here interpreted as SC islands - SCIs), 
separated by regions of vanishingly small order parame- 
ter. The role of temperature is found 16 to be similar to 
that of disorder. 

Recently wc have incorporated thermal phase fluc- 
tuations beyond the BdG theory, 18 and have demon- 
strated the formation of isolated SCIs within the two- 
dimensional disordered negative-U Hubbard model: in- 
creasing magnetic field leads to the loss of SC phase corre- 
lations between different areas in the sample and an even- 
tual SIT. The island picture is very helpful in explain- 
ing some of the experimental puzzles. The coexistence 
and competition between Cooper pair and electron (or 
quasi-particle) transport give rise to the non-universality 
of the critical resistance. It can also explain 19 the non- 
monotonic magneto-resistance behavior observed on the 
insulating side of the SIT. 20 The theory also predicts 
the persistence of the SC islands into the insulating 
phase. Such SC correlations in the insulating phase have 
been observed experimentally by using AC transport 
(and measuring the superfluid stiffness on the insulat- 
ing phase) 21 and by Aharonov-Bohm interference (where 
oscillations with period h/2e in the magnetoresistance 
were observed in the magneto-resistance on the insulat- 
ing phase, indicating the presence of Cooper pairs). 22 

In this picture, the approach to the SIT from the insu- 
lating side is governed by the emergence of phase cor- 
relations between the islands as temperature or mag- 
netic field are lowered, eventually resulting in a macro- 
scopically coherent sample. This picture is in agree- 
ment with the (quantum) percolation scenario. 23 ' 24 Ev- 
idence for percolation-like behavior is substantiated in 
several experiments. 8 ' 9,11 ' 25, 26 Interestingly, the theory 
also demonstrated that phase correlations between dif- 
ferent islands may be explored even within mean-field 
theory, by looking at the response to parallel (Zeeman) 
magnetic field (see below). Thus, in principle, the SIT 
can also be investigated by solving the BdG equation in 
a finite magnetic field. In this work, however, we focus 
on the effects of magnetic field on the formation and evo- 
lution of SC islands in highly disordered two-dimensional 
SC systems in a wide range of temperatures. (We sep- 
arate the effects of magnetic field into orbital field and 
Zeeman field, which wc sometimes refer to as perpendic- 
ular and parallel fields, even though each of the latter 
leads to both orbital and Zeeman effects.) 

The structure of this paper and the main achievements 
are listed below. In section II the model and method 
of calculations are briefly described. In section III A 
the role of perpendicular magnetic fields is studied (as- 
sumed to act on the orbital degrees of freedom as men- 
tioned above). In particular, it is found that, as the 
field is increased, the size of the SC islands shrinks and 
the strength of the local order parameter diminishes. It 
is explicitly shown that vortices are formed and their 
disorder-induced pinning is analyzed. Moreover, these 
vortices tend to accumulate in places where, in the ab- 
sence of magnetic field, the order parameter was low (that 



is in the "normal regions" of the sample). The presence 
of these vortices suppresses the Josephson couplings be- 
tween the SCIs, leading eventually to completely decou- 
pled islands. Section IIIB is devoted to examining the 
role of a parallel magnetic field, manifested by a Zee- 
man splitting. It is found that at low temperatures and 
weak orbital magnetic fields, weakly disordered systems 
exhibit an abrupt vanishing of the SC order parameter 
in response to increasing parallel field, resembling a first 
order transition, as expected for a uniform system. 27 On 
the other hand, when either temperature, disorder, or 
the orbital field are increased, the system undergoes a 
series of transitions, as a function of parallel field, cor- 
responding to the sequential vanishing of A on separate 
islands, in agreement with Rcf. 18. The emergence of 
local magnetic order in areas where superconductivity is 
destroyed is demonstrated. In Section IV we summarize 
and discuss the results. 



II. MODEL AND BASIC EQUATIONS 

In this section, the model and method of calculations 
are briefly introduced. The starting point is the effective 
tight-binding BdG Hamiltonian 17 on a square lattice, 

HBdG = - X! **j ( c !o- c j<r + H - c -) + X^ e « ~ fii)? 1 ™ 

<Lij>,(j i 

+ ^[A(r i )ct T ct i +A*(r i )c iT c ii ] , (1) 

i 

where the sum < ij > runs on nearest neighbors, ti are 
random site energies taken from a uniform distribution 
[— W, W], with W being the disorder strength, fii = fi + 
^■{fii) is the local Hartree-shifted chemical potential, 

A(r i ) = -|E/|( Cii c iT ) (2) 

is the local pairing potential and U < is an effective 
electron-electron on-site attractive interaction. 17 The or- 
bital magnetic field enters via the nearest-neighbor hop- 
ping integral tjj . In the Landau gauge it can be written 
within the Peierls substitution rule as 

Uj = t Oip(-iay.i) (3) 

if r; — Yj is in the x-direction, and tij = to otherwise. 
In Eq. (3) a = 4>/4>o where <j) is the magnetic flux per 
plaquette, 0o = hc/e is the flux quantum, and t/i is the 
y-direction coordinate (in units of lattice spacing) . Here- 
after, to is set to unity and serves as reference to all other 
energy scales. Using a Bogoliubov transformation one 
can diagonalizc H.BdG exactly and obtain the BdG equa- 
tions, 

f i AW U n {Ti) \ _ ( U n {Vi) \ . , 

V A* -t ) \ Mn) )~ \ «n(ri) ) ■ [ ' 

In Eq. (4), ^u„(r,) = - ^ t ili+s u n (ri + 5) + (e» - 
fLi)u n (Yi) where 8 = ±x, ±y, and Au re (r ! ) = A(ri)u„(rj) 
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and similarly for u„(rj), and the energies are the quasi- 
particle (QP) excitation energies E n > 0. The pairing 
potential A(r^) and the electron density per site rii are 
given, in terms of the QP amplitudes u(ri) and v(ri), 
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A(r 4 ) = |[/(r,)|^u„(r 4 )<(r0 

71 

(m) = 2]>>„(r 4 )| 2 • 



(5) 



In what follows, Eq. (4) is solved numerically on a 
finite (albeit large) two-dimensional square lattice with 
hard- wall boundary conditions, starting with some initial 
guesses for A(r,) and rij. A(r^) and n, are then deter- 
mined via Eq. (5), and the procedure is repeated until 
self-consistency is obtained on each site. In order to be 
able to compare the same system at different magnetic 
fields, the chemical potential is altered in every iteration 
such that the total average density per site (n) = i m 
remains constant (N being the total number of sites in 
the lattice). 



III. NUMERICAL RESULTS 



A. Effect of orbital magnetic field 




FIG. 1: Spatial distribution of the order parameter amplitude |A| 
for different values of the orbital magnetic field, 4>/(fro = 0, 0.05, 0.1. 
Bright regions, corresponding to large values of |A|, constitute "su- 
perconducting islands", separated by regions of small |A| (dark re- 
gions). Increasing magnetic field leads to attenuation of the SC 
order parameter and even to the vanishing of SC order in some 
areas of the sample. Calculation is performed on a 12 X 12 size 
lattice with electron density (n) = 0.875, disorder strength W = 1 
and interaction strength U = 2 (this value for U is maintained 
throughout). 



the magnetic field. As can be seen, the distribution func- 
tion shifts from being peaked around some finite initial 
value to being peaked at A = 0. However, the distribu- 
tion tails persist even at high-field, indicating that SC 
islands arc present. 



We start by examining the effect of an orbital mag- 
netic field on the pairing potential A. In Ref. 19 it was 
argued that the SCIs diminish in size and may even van- 
ish as the magnetic field is increased. This assumption is 
indeed confirmed by the numerical results. Figure 1 de- 
picts the spatial distribution of the order parameter am- 
plitude, |A(r)|, in the two-dimensional system, for three 
different values of the magnetic field. At zero magnetic 
field the order parameter fluctuates in space due to dis- 
order, assuming larger values (bright areas in Fig. 1) in 
places where the effective local potential is small, i.e., 
€i — fii « 0, 15 ' 28 which facilitates fluctuations between 
and 2 electrons on the site. These areas are separated by 
regions with small order parameter (dark areas in Fig. 1) 
and hence constitute "superconducting islands" (as will 
be shown in following sections). As the magnetic field 
is increased, the size of these islands and the strength 
of the order parameter diminishes, and some SCIs vir- 
tually disappear (that is, at very strong fields the order 
parameter there becomes smaller than the numerical tol- 
erance). It should be pointed out that while the average 
value of the order parameter diminishes by nearly 50% 
when the flux is increased from 4>/4>o — to 4>/<po = 0.2, 
the maximal value of the order parameter drops by less 
than 10%, meaning that the concept of "SC islands" per- 
sist well within the strong field region (it is the loss of 
correlations between the islands that drives the system 
into its normal state). 

The distribution function of the order parameter am- 
plitude, P(|A|), is plotted in Fig. 2 for different values of 
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FIG. 2: (Color online) Distribution function of the SC order pa- 
rameter amplitude |A|, for seven different values of magnetic field, 
4>/<Po = — 0.18 (calculated for a system size 16 X 16, averaged 
over the entire lattice and over 25 realizations of disorder). With 
increasing magnetic field the distribution shifts from being peaked 
at finite value (the uniform-system value) to being peaked at zero 
value. Note, however, the persistence of tails at large |A| even at 
high magnetic field, indicating the existence of regions with large 
|A|. Inset: the average order parameter amplitude, (|A|), as a 
function of <f>/<f>o, showing a decrease in the SC order parameter 
with the applied orbital magnetic field. 

In the inset of Fig. 2 the average order parameter am- 
plitude (|A|) is plotted as a function of 4>/<t>o • We note 
that (i) in a finite lattice (|A|) never strictly vanishes 
and (ii) for a given disorder realization, (|A|) may be a 
non-monotonic function of 4>/4>q. In fact, in a clean sys- 
tem (|A|) has been shown to oscillate as a function of 
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magnetic field. 29 Disorder, however, tends to smear the 
non-monotonicity, and upon averaging over disorder real- 
izations one obtains a monotonically decreasing function. 

Another illustration of the interplay between disorder 
and magnetic field is manifested in the location of the 
vortices. In a clean system one expects the vortices to 
form a regular Abrikosov lattice. However, strong disor- 
der modifies the position of vortices, as they are pinned 
within regions where the amplitude of the order param- 
eter has a small value. In Figure 3 the phase of the SC 
order parameter is plotted for different values of the mag- 
netic field corresponding to the consecutive insertion of 
one additional flux quantum into the system. As can be 
seen, a new vortex enters the system each time an addi- 
tional flux quanta is inserted, and vortices do not form a 
lattice. Rather, their position is affected by disorder. 



L 2 x(]>/4>0=0 L 2 x<M)0=l 




FIG. 3: Spatial mapping of the SC order parameter phase for 
values of magnetic field corresponding to a different number of flux 
quanta in the system (i.e., 4>/4>o X (L X L) = 0, 1, 2, 3), on a 40 X 40 
size sample. Arrows point at the positions of the vortices. For each 
number of flux quanta, a similar number of vortices is formed, their 
position is irregular and determined by the disorder potential. 



These correlations can be exposed more clearly. In 
Figure 4 the order parameter amplitude at zero field is 
plotted (Fig. 4(a)), together with the phase of the order 
parameter (Fig. 4(b)), for the same system as in Figure 3 
with two flux quanta. As can be seen, the vortex cores 
are positioned where the amplitude is low (dark regions 
in (a)), i.e., they arc "pinned" to regions with low order 
parameter. The emerging interactions between vortices 
still play a significant role, as, for example, the second 



vortex that enters the system changes the position of the 
first one. 

In Ref. 18 it was shown that the SIT is driven by loss 
of phase-coherence between the SC phases on different 
islands. The fact that vortices appear in regions of low 
order parameter amplitude, i.e. in between the SC is- 
lands, is consistent with the picture that orbital magnetic 
field causes loss of phase coherence between the islands: 
the effective Josephson coupling between the islands is 
reduced with increasing magnetic field, until the Joseph- 
son energy becomes smaller than temperature (or the 
amount of quantum phase fluctuations) and the islands 
decouple. 



A=0 A=0.65 




(a) (b) 



FIG. 4: Comparison between (a) the order parameter amplitude 
at zero field and (b) the phase of the order parameter at finite field 
(two flux quanta), indicating that the vortices are formed at the 
position of weak superconducting order (dark regions in (a)). 



B. Evolution of superconducting islands in a 
parallel magnetic field 

In the previous subsection the roles of disorder and 
orbital magnetic field in generating fluctuations in the 
amplitude of the SC order parameter have been demon- 
strated. However, the question remains whether these 
fluctuations may be regarded as separated SC islands. 
Here we demonstrate how the Zeeman effect resulting 
from an application of a parallel magnetic field helps to 
answer this question. This is carried out by extending 
our previous calculations, 18 and include the electron spin 
through the Zeeman effect. 

When a parallel magnetic field is applied, the single 
electron levels split. It was shown a long time ago 27 that 
upon applying a parallel field, a uniform superconduc- 
tor undergoes a first-order transition. The reason for 
this transition is simple: superconductivity requires oc- 
cupation of paired spin up and spin down states, which 
becomes costly with increasing Zeeman energy. For a 
single pair, one expects the transition to take place when 
the gain from superconductivity, Ao, (the bulk value of 
the order parameter), is overcome by Zeeman energy, 
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h c = (]HbB\\. : where B\\ is the strength of the parallel 
field, g being the g-factor for electrons and (i B is the 
Bohr magneton. In fact, when correlation effects are 
taken into account, one finds that the transition from su- 
perconductivity to a magnetically aligned state occurs at 
h c = Ao/\/2- For ultra-small isolated SC grains it was 
shown 30 that a finite size effect (manifested in a level 
spacing of the order of Ao) may smooth the otherwise 
sharp transition. 

In order to account for the effect of parallel field, a Zee- 
man splitting term J2 ia ohc^c^ is added to the Hamilto- 
nian of Eq. (1). The Bogoliubov transformation now has 
to be performed with spin-dependent quasi-particle am- 
plitudes it n ±( r i) and v n ±{Yi) 1 which now obey the spin- 
generalized BdG equations, 31 



0.25 



£±h 

A* 



A 

-i* ± h 



Mu±(r,) 
v„±(r l ) 



(E n ±h) 



u„±(ri) 
v„±(r l ) 



(6) 

Once these equations are solved, the local pairing po- 
tential and occupation are determined by self-consistent 
equations similar to that of Eq. (5). Note that for singlet 
(S-wave) pairing, the different spin-generalized quasi- 
particle amplitudes are decoupled, and the coupling be- 
tween the different spins enters only in the determination 
of the self-consistent Hartrec shift and order parameter. 

For a clean or weakly disordered system and for paral- 
lel field strength h < h c , the eigenvectors are independent 
of h (due to the simple dependence of the energies on the 
Zeeman term) , and hence the solution for the pairing am- 
plitude A(r) and the local occupation is unchanged. On 
the other hand, for h > h c a solution with A(r) = yields 
a lower free energy. Thus, one should observe a step-like 
behavior of the order parameter. In contrast, in a disor- 
dered system, as A(r) varies in space, one also expects h c 
to be space-dependent. For a system composed of sepa- 
rated SCIs, one thus expects that the order parameter on 
each island will vanish when h > h c {i) (where i is an in- 
dex for identifying the corresponding island). Thus, the 
vanishing of the average value of A is expected to exhibit 
multiple steps, each step corresponds to the vanishing of 
superconductivity in a different island. 

In Figure 5 the average order parameter amplitude 
(|A|) is plotted as a function of h for different orbital 
fields. One can clearly see that for small <f> (weak per- 
pendicular magnetic field) , the system undergoes a single 
sharp transition. In contrast, for finite values of <f>, step 
structure of (|A|) is indeed visible. This is seen more 
clearly in the inset of Figure 5 where (| A|) is plotted as a 
function of h at a perpendicular field 4>/4>o = 0.018, with 
much higher resolution in h. At this resolution, the step- 
like behavior is seen even for this high value of the orbital 
field. Thus, large orbital field suppresses SC correlations 
between the islands. They seem to lose SC order sepa- 
rately and successively, depending on the specific value 
of the gap for each SCI. 

To make this visual, the spatial distribution of the pair- 
ing parameter at </>/0o — 0.014 is plotted in Figure 6 for 



|A| 
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0,05 
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FIG. 5: The average order parameter amplitude (| A|) as a function 
of the Zeeman splitting h for different values of orbital magnetic 
field, 4>/4>o = — 0.014. For small 4>/rj>o, the system undergoes 
a first-order transition, which is smoothed as 4>/(fro increases. For 
4>/<fio — 0.004 steps in (|A|) are visible. Inset: the same quantity for 
<t>/<Po = 0.018, but with larger resolution in h. At this resolution, 
steps in the (jA|) are clearly seen. 



h = — 0.012 (Fig. 6(a-h)). The consecutive vanishing 
of the order parameter in different areas with increasing 
parallel magnetic field is clearly visible. It is noticeable 
that even when the order parameter vanishes, as a result 
of the field, in one area of the sample, the other regions re- 
main almost unaffected. The different regions (for which 
the order parameter vanishes at different parallel fields) 
can thus be indeed regarded as separated SCIs. 




FIG. 6: Spatial distribution of the order parameter amplitude |A 
at <f>/(j>o = 0.014 for h = — 0.12 (a-f). Arrows indicate regions in 
which the order parameter vanishes abruptly (SC islands). 



As already mentioned, another indication for the de- 
struction of the SC order may be elucidated by examining 
the local magnetization M(r) = (n^(r) — ri|(r)), which 
develops locally where SCIs are destroyed by the parallel 
field. This is shown in Figure 7, where we plot side by 
side the spatial distributions of the differential variations 
of the magnetization and of the pairing potential, for a 
finite orbital field and for different Zeeman fields (same 
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FIG. 7: Spatial distribution of the differential changes in magneti- 
zation, 6M(r,h) = M(r,h) — M(r,h — 0.03) (columns 1 and 3), and 
in pairing potential amplitude, <5(r,/i) = A(r, h) — A(r, h — 0.03) 
(columns 2 and 4), for different Zeeman fields h = 0.045 — 0.12, 
and for orbital field 4>/4>o = 0.014. The emergence of magnetic 
correlation occurs simultaneously and is highly correlated with the 
vanishing of SC order. 

values as in Fig. 6). It is clearly seen that the magnetiza- 
tion (dark regions in the left panels) develops only where 
SCIs (dark regions in the right panels) shrink. Notice 
that even when the parallel field strength exceeds the 
critical value (h c w 0.24 for these parameters), the mag- 
netization still fluctuates, now due to correlations with 
the local electron density. 



An important and experimentally testable prediction 
of the above calculation is that a strongly disordered SC 
film subject to a parallel field may have both a zero resis- 
tance state and finite magnetization. This is in contrast 
with a clean film, which is either SC or spin polarized. 

IV. SUMMARY AND DISCUSSION 

In this paper the physics related to formation of SCIs 
in disordered two-dimensional SC films and their evolu- 
tion under an applied magnetic field was addressed. SCIs 
arc defined as specific domains exhibiting large value 
of the order parameter amplitude surrounded by other 
domains for which the amplitude of the order parame- 
ter is small. Using the locally self-consistent solution of 
the BdG equations, it was demonstrated that SCIs are 
formed in disordered SC films. Upon an increase of the 
orbital magnetic field, the size of these islands tend to 
decrease, together with the strength of the order param- 



eter. This is accompanied by loss of SC phase correlations 
between different islands. 

In order to verify that SCIs indeed constitute separate 
isolated domains, the application of parallel magnetic 
fields has been included in our numerical simulations. 
For a uniform or weakly disordered sample it results in 
a first order transition into a normal state. In contrast, 
for medium and high disordered systems and in the pres- 
ence of an orbital magnetic field, it is found that upon 
increasing the parallel magnetic field, the average order 
parameter vanishes in a series of steps, indicating that 
each island undergoes a phase transition at its own turn, 
with its own critical parallel field. This substantiates the 
picture of well separated and isolated SCIs. In Ref. 18 we 
have demonstrated a one-to-one correspondence between 
the phase correlations between the islands, at finite or- 
bital field and zero Zeeman field, on one hand, and their 
sequential response to the Zeeman magnetic field, on the 
other hand. This suggests that the SIT can be studied 
within mean-field theory, even though the latter does not 
include phase fluctuations. 

The above results can in principal be addressed exper- 
imentally, by means of local measurements using, e.g., 
scanning tunneling microscopy, which is also sensitive to 
magnetic order. We predict that, as a function of parallel 
magnetic field, the (average) SC order parameter will ex- 
hibit a behavior resembling that of Figure 5, that is the 
transition will be smoothed when the orbital effects are 
dominant. The local SC gap on an isolated island, on the 
other hand, should exhibit a sharp attenuation character- 
istic of a first-order transition. This will be accompanied 
by local formation of magnetic order, so we expect local 
magnetism coexisting with either bulk superconductivity 
(if the SCIs percolate) or local superconductivity (if they 
do not). 

Lastly we point out that some of the phenomena dis- 
cussed here were also observed in high-T c superconduc- 
tors. Local variations in the SC gap have been observed 
by STM measurements. 32 In the underdoped region, it 
was observed 33 that local SC gap persists well into the 
normal phase. Even more relevant is the observation 34 
in these materials of inhomogeneous magnetic response, 
which persists above T c . Thus disorder may play an im- 
portant role in high-T c superconductors, especially in the 
underdoped regime, leading, for example, to the interpre- 
tation of the pseudo-gap, as an average over the spectra 
of SC and non-SC regions. 
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